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Abstract 

o . 

■ Tensor, matrix and quaternion formulations of Dirac-Kahler equa- 
tion for massive and massless fields are considered. The equation 

\ matrices obtained are simple linear combinations of matrix elements 

in the 16-dimensional space. The projection matrix-dyads defining 
all the 16 independent equation solutions are found. A method of 
(T) , computing the traces of 16-dimensional Petiau-Duffin-Kemmer ma- 

trix product is considered. We show that the symmetry group of the 
^ ■ Dirac-Kahler tensor fields for charged particles is 50(4, 2) . The con- 

CN . servation currents corresponding this symmetry are constructed. We 

analyze transformations of the Lorentz group and quaternion fields. 

■ Supersymmetry of the Dirac-Kahler fields with tensor and spinor pa- 
CO . rameters is investigated. We show the possibility of constructing a 
^ cj gauge model of interacting Dirac-Kahler fields where the gauge group 

i ' is the noncompact group under consideration. 

CD . 

1 Introduction 

The important problems of particle physics are the confinement of quarks 
and the chiral symmetry breaking (CSB) [1]. Both problems can not be 
solved within perturbative quantum chromodynamics (QCD). 

One of the promising methods in the infra-red limit of QCD is lattice 
QCD. Lattice QCD takes into account both nonperturbative effects - CSB 
and the confinement of quarks, and provides computational hadronic char- 
acteristics with good accuracy. A natural framework of the lattice fermion 
formulation and some version of Kogut-Suskind fermions [2, 3] are Dirac- 
Kahler fermions [4-11]. The interest in this theory is due to the possibility 
of applying the Dirac-Kahler equation for describing fermion fields with spin 
1/2 on the lattice [4-11]. 
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Recently much attention has been paid to the study of the Dirac-Kahler 
field in the framework of differential forms [11-24]. The author [12] consid- 
ered an equation for inhomogeneous differential forms what is equivalent to 
introducing a set of antisymmetric tensor fields of arbitrary rank. It implies 
the simultaneous consideration of fields with different spins. 

Kahler [12] showed that the Dirac equation for particles with spin 1/2 can 
be constructed from inhomogeneous differential forms. Now such fields are 
called Dirac-Kahler fields. Using the language of differential forms, Dirac- 
Kahler's equation in four dimensional space-time is given by 

(d-5 + m)$ = 0, (1) 

where d denotes the exterior derivative, 5 = — * d* turns n— forms into 
(n — 1)— form; * is the Hodge operator which connects a n— form to a (4 — 
n)— form so, that * 2 = 1, d 2 = 5 2 = 0. The Laplacian is given by 

(d - 5f = - (d5 + 5d) = d^. 

So, the operator (d — 5) is the analog of the Dirac operator •y^d^. The 
inhomogeneous differential form $ can be expanded as 

$ = (p( x ) + ip^(x)dx^ + —ipf, u (x)dx^ A dx"+ 

i ' < 2 > 

-\-—^^p{x)dx 11 A dx v A dx p + —ip^palxjdx^ A dx u A dx p A dx a , 

where A is the exterior product. The form $ includes scalar <p(x), vector 
<Pp{x), and antisymmetric tensor fields if^ u (x), (p^ up (x), (p^ P a(x). The an- 
tisymmetric tensors of the third and forth ranks ip fMup (x), tpp lVpa {x) define a 
pseudovector and pseudoscalar, respectively: 

Ms) = \fT^(?X m = \f vpc v»v P M, (3) 

where e pupa is an antisymmetric tensor Levy-Civita. In fact, the Dirac- 
Kahler equation (1) describes scalar, vector, pseudoscalar and pseudovector 
fields. Some authors (see e. g. [6, 18-21]) showed that Eq. (1) is equivalent 
to four Dirac equations 

( 7 A + m)/ ) W=0, 6=1,2,3,4. (4) 
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The mapping between equations (1) and (4) makes it possible to describe 
fermions with spin 1/2 with the help of Eq. (1), i.e., boson fields ! As we 
have already mentioned, this possibility is used in the lattice formulation of 
QCD and for describing fermions with spin-1/2. 

It should be noted that Ivanenko and Landau [25] considered (in 1928 ) 
an equation for the set of antisymmetric tensor fields which is equivalent to 
the Dirac-Kahler equation (1). Similar equations were discussed after [26-32] 
long before the appearance of [12], and later [33-42]. The author [34, 37] 
found the internal symmetry group 5*0(4,2) (or locally isomorphic group 
SU {2,2)) of the Dirac-Kahler action and the corresponding conserving cur- 
rents. The Lorentz covariance of the Dirac-Kahler equation are also shown. 
The transformations of the 5*0(4, 2) group mix the fields with different spins 
and do not commute with the Lorentz transformations. Later others [36, 39, 
18] also paid attention to this symmetry. The transformations of the Lorentz 
and internal symmetry groups discussed do not commute each other. So, 
parameters of the group are tensors but not scalars as in (the more common) 
gauge theories. This kind of symmetry is also different from supersymmetry 
where group parameters are spinors. The difference is that in our case the al- 
gebra of generators of the symmetry is closed without adding the generators 
of the Poincare group. However the indefinite metric should be introduced 
here. The localization of parameters of the internal symmetry group leads 
to the gauge fields and field interactions. 

The paper is organized as follows. In Section 2 we investigate the ten- 
sor and matrix formulations of the Dirac-Kahler equation for massive and 
massless fields. The tensor and spinor representations of the Lorentz group 
are analized. All independent solutions of equations are found in the form 
of matrix-dyads. It is shown in Section 3 that the internal symmetry group 
of the Dirac-Kahler tensor fields is 50(4,2). For the case of spinor fields 
we come to the ?7(4)-group of symmetry. In Section 4 within the framework 
of the quaternion approach, the six-parameter internal-symmetry subgroup 
and the Lorentz covariance of the Dirac-Kahler equation are considered. 

The quantization of Dirac-Kahler's fields is carried out in Section 5 by 
using an indefinite metric. It is shown in Section 6 that in the field theory 
including Dirac-Kahler fields it is possible to analyze supersymmetry groups 
with tensor and spinor parameters without including coordinate transforma- 
tions at the same time. We show in Section 7 the possibility of constructing 
a gauge model of interacting Dirac-Kahler fields where the gauge group is 
the noncompact group 50(4,2) under consideration. Section 8 contains a 
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conclusion. A method of computing the traces of 16— dimensional Petiau- 
Duffin-Kemmer matrix products is considered in Appendix A. Appendix B 
devoted to the Lorentz transformations and the quaternion algebra. 

We use the system of units h — c — 1, a — e 2 /An = 1/137, e > 0, and 
Euclidean metrics, so that the squared four-vector is v 2 = v 2 + v J = v 2 — v 2 
(v 2 = v\ + v\ + vl v 4 = iv ). 

2 Tensor and matrix formulations of Dirac- 
Kahler equation 

It is easy to show that Dirac-Kahler equation (1) with definitions (2), (3) is 
equivalent to the following tensor equations 

<9„</V - d^cp + m 2 y^ = 0, d u (p^ - + m 2 <^ = 0, (5) 



is the dual tensor, e^ va p is an antisymmetric tensor Levy-Civita; S\2u = 
—%. It should be noted that Eq. (7) is the most general representation for the 
antisymmetric tensor of second rank in accordance with the Hodge theorem 
[43] (see also [44-46]). 

If the ip, ip, ifin, ifn and (p^ u are complex values, Eqs. (5)-(7) describe 
the charged vector fields. These equations are the tensor form of the Dirac- 
Kahler equation (1) which was written in differential forms [12]. 

Now we show that in matrix form equations (5)-(7) can be represented 
as the Dirac-like equation with 16 x 16 -dimensional Dirac matrices. The 
projection matrix-dyads defining all the 16 independent equation solutions 
will be constructed [35]. 

The matrix form facilitates the investigation of the general group of inter- 
nal symmetry. To obtain the matrix form for both the massive and massless 
cases, generalized equations are introduced: 



9^,, -ip = 0, d^f, -(p = 0, 



(6) 
(7) 



where 



1 



(8) 



Q^ix + m 2 ip = 0, 
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- d^y - e^apdaTpp + m 2 VV] = 0, 
<9 M W + m 2 ^ = 0. 



(9) 



With m 1 = m 2 = m, ip = -ip, ip^ = rrup^, Vvi = W> = 
■00 = — = lEfivap (ei234 = 1), we arrive at the Dirac-Kahler equations 
(5)-(7). In the case m\ = 0, m 2 7^ (where m 2 is the dimension parameter), 
Eqs. (9) are the generalized Maxwell equations in the dual-symmetric form 
[47]. Let us introduce the 16— component wave function 

f(x) = {M, A = 0, t i,\ M v],ji,0, (10) 

where ip~ = ip^, ifo = ^o- It is convenient to introduce the matrix e^ 4 '- 8 [48] 
with dimension 16 x 16; its elements consist of zeroes and only one element is 
unity where row A and column B cross. Thus the multiplication and matrix 
elements of these matrices are 



.A,BC,D 



T A,D 



$bc, 



CD 



8ac8bd, 



(11) 



where indexes A,B,C,D = 1,2,..., 16. Using the elements of the entire 
algebra e A ' B , equations (9) take the form 



e H,\iiv\ _|_ £ \pu],n _|_ £ u,0 _|_ £ 0,u _|_ £ i/,0 _|_ 5 0,^_|_ 



^ £ M>M _|_ £ [pw],/i 



+ 



(12) 



+ 



mi + e*") + m 2 ( e°>° + - £ H.M + e °> 



*b(x) = 0. 



Let us introduce the projection matrices 



p — £ v,v _|_ £ /*,/^ p — £ ofi _|_ I £ [H.[H _|_ £ o,o 

2 



(13) 



with the properties PP = PP = 0, P + P = 1^ ; the J 16 is the unit 
16 x 16— matrix and 

T u = e^ u] + £ [H '" + e + e '" + e*>° + + (e^ [puj] + e^-' 1 ) . (14) 
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Then Eq. (12) takes the form of the relativistic wave equation 

(r u d v + m 1 P + m 2 P)^(x) = (15) 

which includes both the massive and massless cases. The 16 x 16— matrix 
T v can be represented in the form 

r„ = = /#> + $ , /?<"> = p® + Pi 0) , 

0(1) = £ ^u] + £ [ H ,m p0) = (f^\ + £ \M^ ) (i 6 ) 

Matrices (3^\ (3^ and f3^\ ffi realize 10— and 5— dimensional irreducible 
representations of the Petiau-Duffin-Kemmer [49-51] algebra 

^(1,0)^1,0)^(1,0) + ^(1,0)^(1,0)^,0) = W (l,0) + ^(1,0) (17) 

and Z^ - ** are 16— dimensional reducible representations of the Petiau- 
Duffin-Kemmer algebra [28-31]. These matrices obey the Petiau-Duffin- 
Kemmer algebra (17) and the matrix T u is a 16 x 16— Dirac matrix with 
the algebra: 

+ r M r, = 2V- (is) 

For the massive case when toi = to 2 = to, Eq. (15) becomes 

(T v d v + to) = 0. (19) 

The 16— component wave equation in the form of the first-order Eq. (19) 
was also studied in [33]. Now we find all independent solutions of Eq. (19) 
in the form of matrix-dyads. In the momentum space Eq. (19) becomes 

- ipty p = em^ p , (20) 

where p = T u p u , and parameter e — ±1 corresponds to two values of 
the energy. From the property of the Dirac matrices, Eq. (18), we find the 
minimal equation for the operator p: 

(ip + to) (ip — m) = 0. (21) 
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According to the general method [52, 53], the projection operator ex- 
tracting the states with definite energy (for particle or antiparticle) is given 
by 

m — i£P ^ > 

= -sr- < 22 > 

This operator has virtually the same form as in the Dirac theory of par- 
ticles with spin 1/2. This is because the algebra of the matrices (18) coin- 
cides with the algebra of the Dirac matrices 7 M . However, here we have the 
wave function ty(x) which is transformed in the tensor representation of the 
Lorentz group. It is also possible to use equation (19) to describe spinor par- 
ticles. In this case the wave function ty(x) will be a spinor representation of 
the Lorentz group and equation (19) is the direct sum of four Dirac equations 
(see (4)). This case is used for fermions on the lattice [4-11]. 

Now we consider the bosonic case. The generators of the Lorentz group 
representation in the 16— dimensional space of the wave functions fy(x) are 
given by (see [28-31, 33]) 

•V = - A (r M r„ - r,r M + ryv - rjV) , (23) 

where the matrices T u also obey the Dirac algebra (18) and have the form 
(see (16)): 

r„ = 4 +) -#->. (24) 

It may be verified that the matrices T M and Y v commute each other, i.e. 

r M ,r,]=o. (25) 

The spin projection operator here is given by 

°V = ; e abcPaJbc = ^abcPa {^b^c + ^b^c) (26) 

2 | p | 4 | p | v ' 

which satisfies the following equation 

a p (a p -l)(a p + l)=0. (27) 

In accordance with [52, 53] the corresponding projection operators are 
given by 

S(±l) = -zCFp (<7 p ± 1) , 5(0) = 1 - <7p. (28) 
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Operators £(±1) correspond to the spin projections s p = ±1 and 5(0) to 
Sp = 0. It is easy to verify that the required commutation relations hold: 
Sf±i) — ^(±1)) S(±i)S(o) — 0) Sfo) — S(o). The squared Pauli-Lubanski vector 
a 2 is given by 

° 2 = (J^ £ ^ a P Pl ' Ja p) = ~2 (j2 J l^ 2 ~ J v«r J ™P»P») ■ ( 29 ) 
It may be verified that this operator obeys the minimal equation 

a 2 (a 2 - 2) = 0, (30) 

so that eigenvalues of the squared spin operator a 2 are s(s + 1) = 
and s(s + 1) = 2. This confirms that the considered fields describe the 
superposition of two spins s = and s — 1. To separate these states we use 
the projection operators 

u 2 u 2 
S(o) = 1 ^"5 ^(1) = ~Y (31) 

having the properties Sf 0) S 2 (1) = 0, (S ( 2 0) ) 2 = Sf 0) , (Sf^f = Sf 1} , Sf 0) + 
S 2 ^ = 1, where 1 = J 16 is the unit matrix in 16— dimensional space. In accor- 
dance with the general properties of the projection operators, the matrices 
S? y S?^ acting on the wave function extract pure states with spin and 
1, respectively. Here there is a doubling of the spin states of fields because 
we have scalar ip , pseudoscalar ip , vector ip^, and pseudovector tp^, fields. 
To separate these states it is necessary to introduce additional projection 
operators. We use the following projection operator 

— m — iep 

M '=^r (32) 

which has the same structure as Eq. (22) but with the matrix p = Y v p v 
and an additional quantum number e — ±1. Following the procedure [52, 53], 
16 independent solutions in the form of projection matrix-dyads are given by 

a 2 m — iep m — iep 1 . . — 
A £l± i,, = - ■ — — • -a p (a p ± 1) = M> £ , ±1 , £ • * £i±1 ,„ 

AS = £ ' ! ^ • ! ^ ' (l " o*) =* ee -%„ (33) 
e,e 2 2m 2m v v > ' ' v ' 
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(0) (, _ m-tep m-iep , _ 2 x _ (0) (0) 



2 J 2m 2m 

where operators A £: ±^, A^ correspond to states with spin 1 and spin 

projections ±1 and 0, respectively and the projection operator extracts 

spin 0. The wave function ty p (^,±1,? or or is the eigenvector of 

the equations 

—ip^/p = em^/ p , —ip^p = em^/ p , 



a p m p = Sp^p, <j 2 ty p = s(s + 1)%, 



(34) 



where the spin projections are s p = ±1,0 and the spin is s = 1,0. The 
Hermitianizing matrix, 77, in 16— dimensional space is given by: 

rj = r 4 r 4 . (35) 

This matrix obeys the equations 

TfTi = -Tim (i = 1, 2, 3), i]T 4 = r 4 77 

which guarantee the existence of a relativistically invariant bilinear form [48, 
53] 

Tj7# = ^+77^, (36) 

where ty p = \I/ + r 4 r4, and is the Hermitian-conjugate wave function. 
In the spinor case, when Eq. (12) is the direct sum of four Dirac equations, 
generators of the Lorentz group in 16— dimensional space are 

jjw = \ {t,y u - r v r„) , (37) 

and the Hermitianizing matrix is 771/2 = r 4 . Using the unitary transfor- 
mation we can find the representation 1^ = J 4 <g> 7 M , where J 4 is 4 x 4— unit 
matrix, 7^ are the Dirac matrices and <E> means direct product. In this basis 
the matrices T M become r = 7^ <g> J 4 . 

It is convenient also to use equations 

(m — iep) (m — iep) = 2ip ( - ± ^ > (ip^ — emj , (38) 
o P = ^ +) = -r^-Ma^^ = = --L-e^^-)^-), (39) 
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where the sign (+) in Eq. (38) corresponds to the equality e — e, and sign 
(— ) to e = —e. With the help of Eqs. (38), (39), the projection operators 
(33) are rewritten as 



W = ^JP {±> (*«*> - em) a*> ± l) = *<$, ■ 



where p^ = p^pj^. Projection matrix-dyads (40) extract solutions 
which are the solutions of the equations 

- ip (+) ¥^ = i ( e + e) m¥+\ (41) 

Z 

- ip { - ] = - (e - e) m^- } , (42) 

Z 

= s p ¥±\ a 2 ^) = s ( a + (43) 

For the bosonic case with £ = e, Eq. (41) describes the superposition 
of vector and pseudoscalar fields, and Eq. (42) with e = —e describes the 
superposition of pseudovector and scalar fields (see [28-31]). 

Let us investigate the case of massless fields; Eq. (15) with mi = 
becomes 

(r„d„ + m 2 P) W(a;) = 0. (44) 

In the momentum space, the field function is the solution to the 
equation 

B^ k = 0, B = ik + m 2 P, (45) 

where k = T^k^, k^ = and the matrix B obeys the minimal equation 

B(B-m 2 ) = 0. (46) 

The projection operator which extracts the solution to Eq. (45) is 

~m 2 -B 

a = (47) 

m 2 
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with the equality a 2 = a required by a projection operator. The corre- 
sponding spin operators are given by 



o k = -f eofcfca/^W. (48) 

We have mentioned that the theory under consideration involves the dou- 
bling of the spin states of particles, because there are vector and pseudovector 
fields. To separate these states in the massless case we can use the following 
projection operators 



A £ = ^(i + e r 5 ), 



where T 5 = r 1 r 2 r 3 r4, e = ±1. The matrix A e commutes with the matrix B 
of Eq. (45) and with spin operators (48), and possesses the required relation 
A 2 = A e . As a result, the projection matrix-dyads, corresponding to the 
generalized Maxwell field after extracting spin and spin projections ±1, 
take the form 

n?> = — (i - {m 2 -B)K = ¥^.^\ 

? V ' (49) 

n( ±1} = —o k K ± 1) (m 2 -B)K = tfW • tff > 
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Let us consider the case of spinor particles when the wave function ty(x) 
realizes the spinor representation of the Lorentz group with the generators 
(37) and Hermitianizing matrix 771/2 = T 4 . In this case the wave function 
ty(x) represents the direct sum of four bispinors and the variables tp , ip^, 
^Vi/]) Vv> which comprise ty(x) (10), are connected with components 
of spinors. Under the Lorentz transformations with generators (37) these 
variables do not transform as tensors. Thus the equations for the eigenvalues 
and the spin operator are: 

-ip*^ = em^' 2 \ a^)^) = 8p *W\ 

(50) 

aJ 1/2) = -— : Te abc p a T b T c , 

4 I p I 

where s p = ±1/2. As there is a degeneracy of states due to the 16 - 
dimensionality we should use additional equations with the corresponding 
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quantum numbers. Taking into account Eq. (25) we can use the following 
additional equations to separate states of spinor fields: 



(51) 

4 1/2) = -r^e a6c p a r b r c 

4 I p I 

with s p = ±1/2. We can treat the additional quantum number s p as the 
"internal spin" because matrices (24) obey the Dirac algebra 

rjv + ry\, = 2v (52) 

Thus it is easy to find all independent solutions of equations (50), (51) in the 
form of matrix-dyads: 



A - - - m ~ l£p . m ' t£p .(± + 2s a^ 2 A ( l - + 21 



(53) 



where ^ £ ,e,s p ,s p = ^t„s,s p ,sp^^ £ = s p = ^V^, and we introduce two 
additional internal quantum numbers e = ±1, s p = ±1/2. In [33] the author 
also used a similar construction for the solutions of the field equations but 
without the dyad representation (53). The dyad representation is essential as 
all quantum electrodynamic calculations can only be done using matrix-dyads 
[52, 53]. The necessary method of computing the traces of 16— dimensional 
matrix products is considered in Appendix A. 

3 The Lorentz covariance and symmetry group 
0(4.2) of charged vector fields 

Let us consider the Lorentz group transformations of coordinates: 

= L^ v x' u , (54) 
where the Lorentz matrix L = {L^ u } obeys the equation 

Lfj, a Lua = b^v (55) 
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Under the Lorentz coordinates transformations (54) the wave function 
(10) transforms as follows 



#V) = TV(x), (56) 

where 16 x 16— matrix T realizes the tensor or spinor representations of 
the Lorentz group. Then the wave equation of the first order (19) is converted 
into 

(r^; + m) *V) = (VpL^dv + m) TV{x) = 0. (57) 

We took into account that at the Lorentz transformations (54) the deriva- 
tives dn become d'^ = L^ u d v . The Lorentz covariance of the Dirac-Kahler 
equation (57) occures if the equation 

T M TV = TT V (58) 

is valid. The infinitesimal Lorentz transformations (54) are given by the 
matrix 

where six parameters define three rotations and boosts. At the same 
time the matrix T at the infinitesimal transformations (59) can be written 

as 

T = Jig + -e^J^, (60) 

where Iiq, J^ u are the unit matrix and generators of the Lorentz group 
in 16— dimensional space, respectively. With the help of equations (59), (60) 
and using the smallness of parameters we arrive from Eq. (58) at (see 
[54]) 

r fxJctv Jau^ \l $afjJ^u ^upJ^a- (^1) 

It is easy to verify that generators (23) for bosonic fields and generations 
(37) for fernionic fields obey Eq. (61). This means that Eq. (19) is covariant 
16— dimensional Dirac-like wave equation which can describe as bosons as 
fermions. For the bosonic case the wave function ty(x) is given by Eq.(10) 
but for the fermion case it is a direct sum of four bispinors (see Eq. (4)). At 
the finite Lorentz transformations the wave function transforms according to 
Eq. (56) with the matrix 

T = exp Qs^J^ . (62) 
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We will show that the group SO(4, 2) is the symmetry group of the Dirac- 
Kahler charged vector fields. This will be obtained using the Dirac matrix 
algebra and the minimality of the electromagnetic interaction [37,36]. 

The interaction with electromagnetic field is introduced by the substi- 
tution <9 M — > = dn — ieA^, where is the vector -potential of the 
electromagnetic field. Consider the Lagrangian 



£ = - 



1 



4* r ' H " 



where 

h 



(63) 

<9 M + ieA^, ^ = ^F^T^, and the arrows above the 
show the direction in which these operators act; = d^A u — d u A^ is 
the strength tensor of the electromagnetic field. From the variation of the 
Lagrangian (63) on wave functions \& and we find equations for Dirac- 
Kahler vector fields in the external electromagnetic fields 



(r^w + m) m(x) = 0, (r M Sw - m) = 0. 



(64) 



From the variation of Eq. (63) on the vector-potential A M we get Maxwell 
equations, in which the source is the electromagnetic current = ie^/T^. 
Let us consider the set of 16 linear independent matrices: 



/ = he, 



i — r 

1 H 1 Hi 



/ = r 5 = r 1 r 2 r 3 r 4 , 



= rsr^. 



(65) 



They commute with the operators (T M £>{-) + m) , (T^D^ - mj of Eqs. 

(64) and generate the algebra of the symmetry of Eqs. (64). The algebra of 
the generators (65) is isomorphic to the Clifford algebra with the commuta- 
tion relations: 



[la/3, Ijxv\ — bpulav + ^avlpn ~ 5/3 u I afl 
[I Hi lap] = fifialfi — b^Iai 



I Hi 1 a/3 — S^alp — 8^/3! a , 



1 Hi ly 



215 



4/, 



-47. 



Hvi 



-21, 



Hi 



I Hi I 



-21 



Hi 



la/3-, 7 



(66) 
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Let us introduce the ant i- Hermit ian generators [55]: 



h — ^16) -^56 — — ^65 — 

_1. _ < 67 > 

With the help of Eqs. (67), the commutation relations (66) take the form 
[Iab, Icd] = $bc!ad + $ad!bc — Sac^bd — SbdIac, 

(68) 

[I AB ,I } = 0, A,B,C,D,= 1,2, ...,6. 

The algebra (68) corresponds with the direct product of the group of 6— 
dimensional rotation 5*0(6) and the unitary group U(l) (for real group pa- 
rameters). This group is isomorphic to the £7(4) group. The transformations 
of the corresponding group are given by 

*'(a:) = UV(x), 

(69) 

U = exp [la + J M /3 M + + I^S^ + l£) , 

where a, (3^, u^ v , <5 M , £ are the group parameters; if these parameters are 
complex, we have the GL(4, c) group. For the neutral Dirac-Kahler fields, 
the transformations (60) should leave real components as real components 
with the conditions: a* = a, (5* m = /3 m , (3% = -/3 4 , uj* mn = uj mn , u* mi = -u mA , 
S^ = —S m , 51 = <5 4 , £* = — £ corresponding to the 50(3, 3)<g>GL(l, i?) group. 
Such a contraction of the GL(A, c) group is a consequence of charged fields 
being described by complex fields having more degrees of freedom. 

The requirement that the Lagrangian (63) is invariant under the trans- 
formations (69) leads to the constraints: a* = —a, fi* m = j3 m , fi\ = —04, 
u* mn = u mn , u* mA = -UmA, S* m = S m , SI = £* = f which corre- 

sponds to contraction of the 50(4,2) <g> U(l) group with 16 parameters 
[36,37]. This occurs only for charged Dirac-Kahler fields. The subgroup 
U(l) is the known group of gauge transformations: fy'(x) = exp (la) ^(x) 
(a* = —a) which gives the conservation law of four-current. For neutral 
fields, the group leaving the Lagrangian invariant under the transformation 
is 50(3, 2) with generators J M , 1^ and corresponding parameters ijj* mn = u mn , 
^4 = _w m4, = ~ S m , SI = S 4 . The generators (see (66)) with param- 
eters LO* mn = oo mn , oj* n4 = —ou mi correspond to the subgroup 50(3, 1). The 
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one-parameter subgroup of the Larmor transformations with the generator 
I was mentioned in [28-31]. Only generators of Larmor and gauge transfor- 
mations commute with the Lorentz group generators (23). This means that 
in the general case, the transformations of the group with internal symmetry 
5*0(4,2) do not commute with the Lorentz transformations. The transfor- 
mations of the Lorentz group realize the operation of internal automorphism 
with respect to the elements of the group considered. As a consequence, the 
parameters of this group are tensors. This is the main difference between 
the group being considered and the usual groups of internal symmetry where 
parameters are scalars. 

As the Lagrangian (63) is invariant under the 5(9(4,2) <g> U(l) group we 
find in accordance with Noether's theorem that the variation of the action is 

5S = j d^xd^ (ty(x)T^(x) - 5^(x)T^(x)) = 0. (70) 

As parameters of transformations (69) are independent we find from (70) 
the conservation tensors: 

j M = w(x)r^(x), = ^(x^^T^ix), = ^(x)r M r 5 r a ^(x), 

_ _ _ _ _ (71) 

= ^(:r)r^(x), e M[a/3] = ^(x)T^ [a T^(x). 

These conservation currents were also constructed in [28-31] without consid- 
eration of the corresponding internal symmetry. Conservation of the currents 
(71) follows from the symmetry-group 50(4, 2) ® £7(1) of the Lagrangian (63) 
(Noether's theorem). For the neutral Dirac-Kahler fields, there is a conser- 
vation of tensors C^ a and Q^\ a 0\ corresponding to the symmetry-subgroup 
50(3,2). Using the matrices r M , F M (16), (24) and wave function (10) V{x), 
ty(x) = ^ + (x)T4T4 i , it is easy to verify that in this case (of neutral fields), 
the currents J M , and R^ a are identically zero. 

For the spinor case with the generators (37), the Lagrangian (63) with 
the conjugate function ^f(x) = ^ + (x)T 4 is invariant under the U(A) - group 
transformation (69) with the parameter constraints: a* = —a, (3* = —(3^, 
= 00 ^ i = £* = ~~ £ • m this case the transformation (69) commutes 
with the Lorentz transformations (see (37). The existence of the additional 
quantum numbers s p and e is connected here with the presence of the group 
577(4). The subgroup U(l) with the parameter a is the well known group 
of gauge transformations giving the conservation of the electric current J M = 
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4 Quaternion form of Dirac-Kahler's fields and 
symmetry 

Dirac-Kahler equations (5)- (7) may be written in quaternion form [56, 57]. 
Within the framework of the quaternion approach, the six-parameter internal- 
symmetry subgroup of the Dirac-Kahler equation is considered. It is shown 
that the Lorentz group is the automorphism group of the group under consid- 
eration. The possibility is investigated to reproduce the potentials and field 
transformations induced by the coordinate transformation in the complex 
space-time [58]. 

Let us introduce the following biquaternions (see Appendix B) 
V = e^, F = F M e M , G = G M e M , 
F m = H m - iE m , F 4 = -ip- itp, G M = (y? M + itp^) , (72) 

Using the algebra of quaternions (see Appendix B) we represent Eqs. 
(5)-(7)as 

VF + m 2 G = 0, F — -VG, (73) 

where V = e M <9^; e M = (e^, — e m ) are the conjugated quaternion elements. 
The Lagrangian (63) with the help of the basis elements of the quaternion 
algebra takes the form 

C = -\ (FF + m 2 GG + F*F* + m 2 G*G*) . (74) 

Equations (73) preserve their form under the following transformations 
of the field variables: 

G -> G' = GD, F -> F' = FD. (75) 

In the general case, transformations (75) define the SL(2, c) group as the 
quaternion algebra isomorphic to the algebra of Pauli's matrices (see Ap- 
pendix B). The Lagrangian (74) is invariant under the transformations (75) 
if the biquaternion D satisfies the equality DD = 1. This condition defines 
a 6— parameter group of internal symmetry 5*0(3,1), which is a subgroup 
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of 5*0(4,2) investigated in the previous Section. We can use the following 
parametrization of the biquaternion D: 

D = exp n, (76) 

where n is the vector biquaternion with six parameters. The complex 
quaternion (76) obeys the equation DD = 1. The finite transformations 
(75) with the biquaternion D (76) define subgroups SU (1,1) and SU(2) 
for the real and complex parameter rii, respectively. It is easy to verify 
that transformations (75) corresponds to (69) at a = /3 M = <5 M = £ = 0, 

W[m4] = -W[4 m ] = Iui 7l m , UJ[ mn ] = e mn feRe Ufc. 

Under the transformations of the Lorentz group, the quaternions G, F, 
V and V are changed as follows 

G L = L*GL, F l = LFL, 

(77) 

V L = L*VL, V = LNL* 

with the condition LL = LL = 1. It is obvious that the field equations 
(73) are invariant under the Lorentz transformations (77) which guarantee 
the relativistic invariance. Now let us consider the transformations of the 
Lorentz group L and D— transformations (75). We have 

G DL = L*GDL, G ld = TGLD, 

(78) 

F DL = LFDL, F LD = LFLD. 

It follows from Eqs. (78) that the Lorentz group does not commute with 
the group of the internal symmetry (75) as G DL ^ G LD , F DL ^ F LD . Under 
the Lorentz group transformations, the biquaternion D transforms as 

D L = LDL. (79) 

It is seen from Eq. (79) that parameters of the group (75) are trans- 
formed under the tensor representation of the Lorentz group. Eq. (79) also 
denotes that the Lorentz group is the automorphism group of the group under 
consideration. 

Let us consider the complex four- dimensional space where coordinates are 
complex potentials G^ = ip^ + iCp^. The group of transformations of four- 
dimensional rotations in this space is homomorphic to the group SO (A, c); 
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it leaves the quadratic form G 2 ^ invariant. We suppose that the space-time 
coordinates are not transformed here. The transformations of this group 
are given by (Appendix B) 

G -> G' = SGR, (80) 

where 55 = RR = 1. Because coordinates are not changed under this 
transformation, the quaternions V and V are also not changed. Transforma- 
tion (80) leaves Eqs. (73) invariant. Indeed, 

F' = -VSGR, 

VF' + m 2 G' = V (-VSGR + m 2 SGR) = 5 (VF + m 2 G) R 
This is extracted three subgroups SL(2,c) from 5*0(4,0): 

5 = 1, RR=1, G' = GR, F' = FR, 

R=l, 55 = 1, G' = SG, F' = -V5G, 

S = if, RR = 1, G' = R*GR, F' = -VR*GR 

It is obvious that the group (82) coincides with the group (75). 

Now we discuss the possibility of inducing the transformations (75) by 
the Lorentz transformations in the complex space-time (see Appendix B). 
Transformations of the complex coordinates z 1 = LzR of the complex Lorentz 
group 50(4, c) induce the following transformations of biquaternions V, V, 

G and F: 

V = LVR, V' = R7L, 

(85) 

G' = LGR, F' = RFR, 

with the constraints: RR = 1, LL =1. At the transition to the complex 
space-time, transformations of the complex Lorentz group 50(4, c) and in- 
duced transformations of biquaternions V, V, G and F Eq. (85) retain 
the invariant form of Eqs. (73). In the case of the ordinary Lorentz group 
50(3, 1), we should set L = R . The transformations of the group of the 
internal symmetry of potentials G, Eq.(75), and the transformations of G, 
Eq. (85), from the complex Lorentz group at L — 1 have the same form; i. e. 
they are not different. But transformations of the field biquaternion F in Eq. 
(75) and (85) are different. Therefore the transformations of the potentials 



(81) 

(82) 

(83) 
(84) 
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iPn , (fn can be induced by the transformations of the complex Lorentz group 
5*0(4,0) but the fields ip^, ip, <p> cannot. 

The possibility of considering transformations (75) for the field variables 
is based on the fact that F 4 = —p — ip 7^ but transformations (85) are valid 
also for F 4 = 0, i. e. for equations without additional scalar and pseudoscalar 
fields. At ip — <p — 0, (ftp — Eqs. (73) represent the quaternion form of the 
Proca equations. For them the transformations (75) are not possible. 

It should be noted that the transition to complex space-time is used in 
the investigation of some general problems of quantum field theory; e. g. the 
solution of some specific problems in electrodynamics [59-60]. 

5 Quantization of fields 

The quantization of Dirac-Kahler's fields will be carried out by using an 
indefinite metric. It will be shown that the renormalization procedure is 
carried out in the same manner as in quantum electrodynamics [61, 62]. 

The Lagrangian of charged Dirac-Kahler fields (63) within four-divergences 
can be written (when electromagnetic fields are absent) as 

C = -H(x) (T^ + m) *(x), (86) 

where ty(x) = ^(rr) 4 "]^!^ corresponds to the tensor representation of the 
Lorentz group, where the 16-component wave function ^ describes scalar, 
pseudoscalar, vector and pseudovector fields. In the case of spinor repre- 
sentation of the Lorentz group, \I/ is the direct sum of four Dirac bispinors, 
fy(x) = ^/(x) + T 4: , and the quantizing procedure is similar to the Dirac theory. 

Now we will consider the case of the boson fields. Using the canonical 
quantization, one arrives at the commutators 

where M, N = 1,2,.. .,16. It follows from (87) that it is necessary to in- 
troduce the indefinite metric, as for finite-dimensional equations (with the 
diagonal matrix T 4 ) which describe fields with integer spins; only fields obey- 
ing the Petiau-Duffin-Kemmer equation have positive energy [63]. With the 
help of Eqs. (10), (14) we get the following commutation relations for the 
tensor fields: 

[ip(x), mip* (x')] t=t , =iS(jc- x'), [<p{x),mp*{x')] t=t , = -i£(x - x'), 



= „ = (r 4 W(x-x'), (87) 
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[™*Pk(x), <P* mn (x% =tl = ie kmn 5(x - x'), 



(88) 



[rmp k (x), <P* n4 (x% =t , = S kn S(x - x'), 

plus complex conjugated relations. In the momentum space, the equation 
of motion of fields with spins and 1, takes the form 

(m±ip)^ ± {p) = 0, (89) 

where p = p^T^, ^(p) are positive (^(p)) and negative (^"(p)) fre- 
quency parts of the wave function corresponding positive p > and negative 
p < energies of particles, respectively. For each value of the energy, there 
are eight solutions with definite spin, spin projection and addition quantum 
number. Wave functions \I/ ± (p), \I/ (p) can be expanded in spin states as 
follows 

*£(P) = «?(pK^ (P), *£(p) = a^(p)v^(p) (r,T 4 ) NM , (90) 

where index s = 0, m,n, (m, n = 1,2,3), and operator a* + (p) is the 
creation operator of a particle in the scalar state (s = 0), pseudoscalar state 
(s = 0), vector state (s = m), pseudovector state (s = n), and aj(p) is 
the annihilation operator of a particle. The normalization conditions for 
solutions (90) are different from the Dirac bispinors case, and are given by 

VN'Hp)(^) NM v r M T (p) = ±e s 5 sr , (91) 
^(pK'^p) = e s 5 sr -, = ^(pjrA) , (92) 

t;*^(p)r 4 ^ ± (-p) = 0, (^(p))* = v* s *(p), (93) 

and e s = 1 at s = 0, m, and e s = — 1 at s = 0, n. We use here the 
normalization on the charge and in the right hand sides of Eqs. (91), (92) 
there is no summation on indexes s . The summation formula on indexes s 
corresponding to the normalization conditions, has the form 

E^ ± (P)^ T (P)=( ! ^) • (94) 

V Z P0 I MN 

If we take the trace of the matrix (94) and compare it with the expression 
(92) summed over all states s, we will get the equality. Multiplying Eq. (94) 
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into the matrix r 4 , and then calculating the trace of both sides of the matrix 
equality, we arrive at Eq. (91). 

The appearance of the coefficient e s = ±1 in the right hand side of Eq. 
(92) reflects the fact that the energy of vector and pseudoscalar states is 
positive, and the energy of pseudovector and scalar states is negative. We can 
also come to this conclusion using the expression for the energy-momentum 
tensor 

T\xv = -^(x)T^(x) (95) 



which follows from the Lagrangian (86). Taking into account the expan- 
sions for wave functions: 



d 3 p 



*(x) = (2tt)- 3 / 2 | [^+(p)e ip * + ^-(p)e- 

= (2vr)- 3/2 | \^ + (p)e ipx + V~(p)e~ 
found from Eq. (95), the energy- momentum- vector 

ty(x)T 4 d^{x)d 3 x, 

and Eqs. (90)- (93), we arrive at the following expression 

P, = /pmE £ . (< + (p)«;(p) + a*-(p)af(p)) d 3 p. 



(96) 



(97) 



(98) 



To have the positive energy in accordance with the Eq. (98), we should use 
the following commutation relations for creation and annihilation operators: 



«r(p),a r + (p') 



a:(p),K + (p') 



e s 5 sr 5(p - p'). 



(99) 



With the help of Eqs. (90), (94) and (99) we find 

*mW,*1(!/)] = "(2-r 3 / {m+ f )MN e-^-y^p 



dx, 



m 



2 P o 
A_(x-y), 



(100) 



MN 



^(x),W N (y)] = (2vr)- 3 / {m ^MN e i P ( x - y)d s p 

J 2p Q 



(101) 
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-y)- 

It follows from Eqs. (100), (101) (by taking into account the definition 
of the Pauli- Jordan function [54, 64] A = % (A + (x) — A_(x))) that 

*(x),*(y)] = S(x - y), S(x - y) = S + (x -y) + S-(x-y), 

(102) 




where the function S(x — y) satisfies the following equations 

(iV A + m^j S(x - y ) = i^- rrfj A (x - y) = 0. (103) 

From Eqs. (102), at t — t', and using Eqs. (100), (101) we arrive at the 
commutator (87). With the help of the relationships (102), we can find the 
chronological pairing of the operators: 

(T* M (x)* N (y)) = S c MN (x-y) 

= Q(x - y )S% IN (x -y)- Q(y - x )S^ N (x - y) (104) 

= % J ( l P- m )MN c ip(x~y) d A p 

(27r) 4 J p 2 + m 2 — ie 

which has formally the same form as in quantum electrodynamics (QED). 
Here Q(x) is the well known theta- function [64] and x is the time. 

It is seen from Eq. (104), that the Feynman rules for particles with spins 
0, 1 interacting with the electromagnetic field, eventually are the same as in 
QED. We should not, however, here use the QED factor 77 = (—1)', where 
I is the number of loops in the diagram due to different statistics [64] . The 
difference is in the number of spin states of the charged particle, and in the 
dimension of matrices T^. As the propagator (104) formally coincides with 
the electron propagator of QED, so all divergences can be cancelled by the 
standard procedure, i. e., we have here a renormalizable theory. All matrix 
elements of quantum processes describing the interaction of particles with 
multispin 0, 1 coincide eventually with the corresponding elements in QED. 
The difference is in the density matrix \1> • \1> which we found in Sec. 2. 

The commutation relations (99) with sign (— ) in the right hand side 
(at e s = —1) require the introduction of the indefinite metric. The space 
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of states is divided into two substates: H p and H n with positive (H p ) and 
negative (H n ) square norm. The vector and pseudoscalar states correspond 
to a positive square norm, and pseudovector and scalar states — to a negative 
square norm. The total space is the direct sum of the two subspaces H p and 
H n . 

6 Super symmetry of Dirac-Kahler's fields 

For the Dirac-Kahler fields, it will be shown that in field theory it is possible 
to analyze transformation groups with tensor and spinor parameters without 
including coordinate transformations at the same time [65]. 

A graduated Lie algebra must be absolutely related to transformations 
involving space-time coordinates [66], but it seems perfectly obvious that this 
is always the case if we are dealing with transformations whose generators 
are of neither a tensor nor spinor nature (see [67, 68], for example). 

A theory of Dirac-Kahler's fields, however, raises the possibility of con- 
structing a transformation group with tensor and spinor generators which 
does not at the same time include coordinate transformations. 

Let us consider the field equations 

+ \ ( m i + m 2^j G(x) + i (m 2 - mi) 7 5 G(:r)75 = 0, (105) 
where 7 M are the Dirac matrices, and the matrix G(x) is 
G(x) = ip (x)I 4 - Vv(z)7a. + ^l^]( x hlfilu] + «Vv(z)7m75 ~ #o(a075- (106) 

The quantities i[)o(x), V'mO 1 ')' i } \}iv\{ x )i ' ) /v( x )> anc ^ $o( x ) i n Eq. (106) are, 
respectively, a scalar, a vector, an antisymmetric tensor, a pseudovector and 
a pseudoscalar; under the Lorentz group, G(x) transforms as follows: 

G L (x) = SG(x)S-\ S = exp Qe^T^]) - ( 107 ) 

where e^ u are the Lorentz group parameters. Multiplying Eq. (105) 
successively by the Clifford-algebra elements 7^4: 2/4, 7^, (1/2)7^7^, 7^75 
and 75, and taking the trace, we find the tensor equations which coincide 
with Eqs. (9) including the massive and massless cases. The case of a 
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massless field corresponds to the choice m\ — 0, while a massive field (at 
mi = m 2 = m ) is described by an equation of the type 

( 7 A + m)G(x)=0. (108) 

The matrix equation (108) is equivalent to the massive Dirac-Kahler 
equation. It is obvious that Eq. (108) describes spinor particles when the 
4 x 4— matrix G(x) represents four bispinors, and that it describes scalar, vec- 
tor, antisymmetric tensor, pseudovector and pseudoscalar fields when G(x) is 
expanded by (106). For the spinor case, however, G(x) transforms as follows 
(in this case we use the index 1/2): 

G 1/2 (x) -> G{ /2 {x) = SG 1/2 (x), S = exp Q^l^]) • (109) 

The Lagrangian corresponding to Eq. (108) is 

C = - l -ti [Gixh^Gix) - Gix^X ^) + 2mG{x)G{x)] , (110) 

where G(x) = j^G^x)^, the arrow specifies the direction in which the 
differential operator acts and tr means the trace of matrices. After taking 
the trace in Eq. (110), we arrive at the Lagrangian which is equivalent to 
Eq. (86). 

Equation (108) is invariant under the following transformations of the 
matrix quantity G(x): 

G(x) -> G\x) = G{x)D. (Ill) 

The relativistic invariance of Eq. (108) is retained if the matrix D trans- 
forms under the Lorentz group as 

D{x) -> D L (x) = SDS~\ (112) 

i. e., if the generators of the transformation group (111) are of a tensor 
nature with respect to the Lorentz group. In the spinor case, all parameters 
of transformation (111) are scalars. Therefore the Lorentz transformations 
(109) and the transformations of the internal symmetry (111) commute each 
other. If we make the Lorentz (109) and symmetry (111) transformations for 
the case of spin 1/2, one gets 

G 1/2 (x) - G? /2 (x) = SG 1/2 (x)D. (113) 
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The commutation of these two groups of transformations is obvious from 
Eq. (113). If we put D = S^ 1 in Eq. (113) we arrive at the law of the 
transformation of the tensor fields (107). That is why it is possible to describe 
the spinor particles with spin 1/2 by the tensor fields using the expansion 
(106). The same conclusion follows from the formalism of Sec. 2. 

The requirement that the Lagrangian (110) be invariant under transfor- 
mations (111) leads to the condition 

DD=1, 7J(a;) = 74^+74. (114) 

Writing D in the form 

D = exp (ial 4 + + ^l^l^lv] + <^7m7s + Pl^j , (U5) 

where U4, 7^, ^[nJv], 7^75 and 75 are the generators of the group GL(A, C), 
we find from condition (114) a restriction on the parameters: a* = a, 

firn = fimi @4 = ~ 0Ai ^Knn = ^mnj ^m4 = _ ^m4i $m = S m , 5^ = —S4, p* = p, 

in accordance with a singling out of the SO (A, 2)®U (1) subgroup (or locally 
isomorphic to [7(2,2)). The subgroup [7(1) corresponds to the gauge trans- 
formations that conserve electric current. So the symmetry group of equation 
(108) is GL(A,C) and the corresponding Lagrangian — 5*0(4,2) <g> [7(1) for 
a case of tensor fields. In the case of spinor fields, the symmetry group is 
[7(4), in accordance with conclusions of Sec. 3. 

From the invariance of the Lagrangian (110) under transformations (111) 
((114) is taken into account) we find conservation laws for quantities of the 
type 

6 M = ^tr (i,G(x) lA G(x) - 7^(20747+74^)) , (116) 

where j A = 7 M , (1/2)7^7^, 7^75, 75 and 7^ are the complex conju- 
gated Clifford-algebra elements. The conservation currents (116) were found 
in Sec. 3 in another formalism. From the physical standpoint, the appear- 
ance of this symmetry results from a mass degeneracy of the spin states of 
the particle which are mixed by transformations (111). The symmetry is 
preserved in a nonlinear generalization of equation (108) (equations of the 
Heisenberg type): 

(7A + m ) G{x) + lG{x)G{x)G{x) = 0, (117) 
where / is the coupling constant. 



26 



The matrix G(x) corresponds to a second-rank bispinor G&. If the bispinor 
satisfies the Dirac equations with respect to each index simultaneously, we 
find a system of Bargmann-Wigner equations [69] which describe a particle 
with spin and system of particles with spin 1. 

By jointly analyzing Eqs. (108) and the Dirac equation 

(tA + "0*0*0 = o ( 118 ) 

which is invariant under phase transformations of wave function ty(x) 
[if(x) — > ty'(x) = exp(i9)ty(x)], we can construct a symmetry group which 
incorporates these phase transformations and transformation (111) as a sub- 
group. The system (108), (118) is invariant under the following transforma- 
tions: 

G\x) =G(x)D + V(x)-(, 

(119) 

V'(x) = *(x)A + G(x)£, 
where *&(x) ■ £ = (^ a (x) ■ C 13 ) is the matrix-dyad, and ( , £ are bispinor- 



parameters, and A is the complex number-parameter. Transformations (119) 
form a group with the following parameter composition law: 

D" = D'D+e-c, c" = a'c+c'a 

(120) 

A" = A'A + C'£, £" = £'A + 

Under the Lorentz group, the parameters ( and £ transform as bispinors 
^ and respectively and are constant quantities, independent of the space- 
time coordinates. In order to preserve the relationship between the spin and 
the statistics, we must require that the parameters £ and ( anticommute: 
— {C^C } = 0. The need for this condition can be seen directly 
from the commutation relations for boson and fermion fields and from the 
explicit form of transformations (119). 

To establish the group structure of transformations (119) it is convenient 
to use a 20— component column function <&(x) whose first components are 
formed by the elements of the lines of the matrix G^(x) (in the order of an 
alternation of lines and of the elements in them). The other four components 
correspond to the wave function ^(x), so 

*<*) = ( S! ) ■ (i2i) 
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A direct check confirms the following form for writing transformations 
(119): 

&{x) = (h®B)*{x), B =[f V)' (122) 
where ( j is a column, and £ is a row and m T ) = I)/*,. Under 

V / a(3 

the condition tr(.D T ) = A, the transformations in (122) correspond to the 
graduated Lie algebra SL(4 | 1), in the notation of [70]. The form in (122) 
for the transformations (119) is of a standard type [71], where D T and A are 
even elements of a Grassmann algebra, and £ and ( are odd elements (i.e., 
anticommuting elements). 

A fundamental distinction between this symmetry and the "ordinary" 
super symmetry, with tensor and spinor generators, is that the corresponding 
superalgebra is closed without appealing to the generators of a Poincare 
group. 

By analogy with the description of fields with a maximum spin of 1, fields 
with a maximum spin of s/2 can be described by the equations 

+ m) G aia2 ... as = 0. (123) 

Dirac equation (118) and Eq. (108) constitute a particular case of Eq. 
(123), with {G Ql ) = ty(x) and (G aiCt2 ) = G(x). A particle with a maximum 
spin of 3/2 and a rest mass is described, for example, by the function G aiCt2a:i . 
These internal-symmetry transformations can be generalized immediately to 
the case of particles with maximum spins of 1 and 3/2, which are the particles 
of the greatest physical interest. 

What possible physical applications could this symmetry have? In an 
analysis of systems consisting of two (mesons) and three (baryons) quarks 
in the s state, and in the absence of a spin-spin interaction, the hadrons 
may be described as multispin particles with a rest mass. In this case the 
symmetry under consideration holds rigorously, and the hadron interactions 
can be associated with the internal-symmetry group 5*0(3, 1) ®SU (3), where 
5*0(3, 1) is the internal symmetry group, which forms with the Poincare 
group a semidirect product. As has been mentioned in the literature [72] 
the latter property is a necessary feature of strong-interaction symmetry 
groups with incorporated quark spin. From this standpoint the considered 
supersymmetry corresponds to an internal symmetry of a field theory which 
incorporates, along with composite particles, their structural components. 
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A further study of this symmetry will be required, of course, to take into 
account its possible extension to interacting fields. 

7 Non-Abelian tensor gauge theory 

We will show the possibility of constructing a gauge model of interacting 
Dirac-Kahler fields where the gauge group is the noncompact group 5*0(4, 2) 
under consideration [73]. 

The starting point in the introduction of Yang-Mills fields is the localiza- 
tion of parameters of the symmetry group, the transformations of which do 
not affect the space-time coordinates. We consider fields ty(x) that possess 
certain transformation properties under the Lorentz group and that may be 
transformed under a certain representation of the internal symmetry group 
(usually compact and semisimple of type SU(n)). For example, QCD consid- 
ers the fermionic fields of spin 1/2, which are transformed under the funda- 
mental representation of SU (3), in which the colored quarks are the principal 
objects. In other words the concept that there exist some internal quantum 
numbers (isospin, colour, etc.) is a requisite physical element of non-Abelian 
gauge theory. However, it has been a long-standing problem to describe par- 
ticles without involving internal ("isotopic") spaces (for example, the works 
of Heisenberg, Diirr [74, 75] and Budini [76]), but using instead an adequate 
generalization of the known relativistic wave equations (RWE) (see reference 
[77] and other references below). It becomes imperative today, for an ap- 
proach of this kind, to imply the concept of a non-Abelian gauge field— the 
carrier of interactions. The possibility of constructing the gauge theory is 
inherent in the theory of RWE. The theory might be based on the concept 
of a multispin or, equivalently, of a particle having several spin states. The 
equation remains coupled— i. e., it describes, as the ordinary Dirac equation 
does, a particle (and antiparticle) having a set of states, rather than a set 
of particles, as is the case with, for example, the equation corresponding 
to the direct sum of the Dirac equations. Transformations of internal sym- 
metry groups result in a mixture of states related to different values of the 
spin squared operator. Their localization leads to non-Abelian gauge fields 
having multispin 0, 1, 2. In this case the dynamics of the interaction of parti- 
cles with multispin 0, 1 are associated with the change of their state through 
the exchange of particles with maximum spin 2. The theory constructed in 
this manner represents a space-time analogue of gauge theory with internal 
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symmetry. 

An attractive possibility is to describe quarks by the Dirac-Kahler equa- 
tions. The authors [4-11] used lattice version of the Dirac-Kahler equations 
describing fermions by inhomogeneous differential forms. This is equivalent 
to introducing a set of antisymmetric tensor fields of arbitrary rank for de- 
scribing the fermion matter fields. As shown in the Introduction, we arrive 
at the Dirac-Kahler formulation which includes a scalar, a vector, an anti- 
symmetric tensor, a pseudovector and a pseudoscalar field. Here we consider 
the continuum case of the equations and introduce non-Abelian tensor gauge 
fields (gluon fields) for interacting quarks. In this point of view, quarks pos- 
sess the multispin 0,1. As we have already shown, in the continuum case 
the equation for the 16— component Dirac equation can be reduced to four 
independent 4— component Dirac equations. We proceed to use the language 
of tensor fields to formulate non-Abelian tensor gauge theory. 

The requirement that the Lagrangian (63) be invariant under local trans- 
formations (69) leads to the necessity of introducing a compensating field 
A*}, where the index B is "internal" (in our case it represents a set of tensor 
indices specifying a scalar, a four- vector, a skew second-rank tensor, an axial 
four- vector and a pseudoscalar). The gauge invariant Lagrangian has the 
known form: 

C = [r„ (d, - ieA, - gAp B ) + m] *(x) - -J% - - (F*f , (124) 



where 



•F l UV ■ 
pB _ o A B 



d^A v 



dvA^, 



9uA^ — -gcQ D A^A^ 



(125) 



and Cq D are the structure constants of 5*0(4,2) group; J 7 ^, F^ u are the 
strengths of the electromagnetic and "gluon" fields, respectively; g is the 
"gluon" coupling constant and B = {fx, [a(3], ji, 0} . The localization of the 
U(l) group produced the electromagnetic field with four-potential A^. 

The corresponding wave equations which follow from the Lagrangian 
(124) are 

d v T^ = J M , (126) 
duF% + 9C B DC A C V F» = J* (127) 

(128) 



r M (d„ - ieA^ - gA B I B ) + m] = 0, 
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where J, = ieVF^, J* = g^T^; I B = r„, (1/4)1^, T 5 , r 5 r M 
(see Eq. (66)). The conservation current for the non-Abelian fields is 

J* = J*-gc B DC A C ;F°. (129) 

In the general case the gauge field multiplets include a second-rank 
tensors A", A", a third-rank tensor antisymmetric over two indices A^ v \ and 
an axial four- vector A®. The structure constants Cq D transform in the tensor 
representation of the Lorentz group. The gauge fields carry the maximal 
spin 2. Indeed, the second-rank tensor A" is transformed on the following 
superposition of the irreducible representations of the Lorentz group: 

G' D ® & I) = (0, o) ® (0, 1} (1, o) (1, i} (i3o) 

corresponding to the fields with the spins 0, 1, 2. The third-rank tensor 
A^ realizes the representations 

(^[(o.^o>H^(^^MM) <-) 

which also contain fields with spins 0, 1, 2; hence we come to the same 
conclusion about the spin of the gauge fields. 

We can also consider the localization of some subgroups of the total group 
5*0(4, 2). Then the Lagrangian (124) will be invariant under the local trans- 
formations of this subgroup. The main requirement is to extract the sub- 
group in a relativistic manner. We suggest the following subgroups and their 
corresponding generators 

50(3, 1) - {V]} , SO(3,2) - {V],/ a } , 

(132) 

SO(4,l)-{/ [H ,/ Q }, GL(l,R)-{l}. 

The possibility of constructing a dynamic theory with all the main prop- 
erties of gauge theories, but based upon notions of space-time rather than 
on new internal quantum numbers, is of evident interest. 

The absolute group symmetry G corresponds to the semidirect multipli- 
cation of the Poincare group P on the internal symmetry group D (5*0(4, 2) 
group): 

G = P-D, (133) 
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and the transformations of the symmetry group (D— group) commute 
with the transformations of the subgroup of four-translations T 4 . Follow- 
ing [78], it is possible to define the "auxiliary" Poincare group P', which is 
isomorphic to P, by the relationships: 

p' = {l;,}-t 4 , l;, = v-v, (134) 

where 1^ are the generators of the internal symmetry group 5*0(3, 1) (see 
(66), (69)). As a result we have 

G = P ■ D = P' <g> D, (135) 

i. e., the absolute group symmetry G represents the direct product of 
the auxiliary Poincare group P' and the internal symmetry group D , taking 

p4, I^v] = 0. In our case of the Dirac-Kahler 



into account that 



T I T 



equation, 1^ = (1/4) (1^1^ - r^) and 

J> = \ (r,r, - r„r M ) . (136) 

The generators of the "auxiliary" Lorentz group, U , commute with the 
generators of the internal symmetry group, Iab (67), and the wave func- 
tion transforms in the spinor representation of the group {L^„}. This con- 
firms that the Dirac-Kahler equation describing a set of antisymmetric tensor 
fields by the inhomogeneous differential forms, can describe spin 1/2 parti- 
cles (pseudoscalar and pseudovector fields are equivalent to the antisymmet- 
ric tensors of fourth and third ranks, respectively). The non-Abelian gauge 
theory under consideration is an analogy to the ordinary non-Abelian gauge 
theory of spin 1/2 particles interacting via gluon fields with internal symme- 
try group 50(6) (or SU (4)). However in our case we have noncompact gauge 
group 50(4, 2) (or SU(2, 2)) which requires the introduction of an indefinite 
metric. 



8 Conclusion 

We have considered Dirac-Kahler equations which can be represented as the 
direct sum of four Dirac equations. The main feature of such scheme is 
the presence of the additional symmetry associated with non-compact group 
in the Minkowski space. The transformations of this group mix fields with 
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different spins and they do not commute with the Lorentz transformations. 
As a result, the group parameters realize tensor representations of the Lorentz 
group. This kind of symmetry differs from the colour and flavour symmetries 
of QCD and supersymmetry. At the same time it was shown that the Dirac- 
Kahler fields allow us to introduce graduated groups with tensor and spinor 
parameters without including coordinate transformations. 

The field scheme considered allow us to construct gauge theories with 
different non-compact groups 0(3,1), 0(4,2), 0(3,3), where "gluon" fields 
carry spins 0, 1, 2. Some of these sup-groups become compact groups in 
the Euclidean space-time. The theory constructed represents a space-time 
analogue of gauge theory with internal symmetry but there is a difficulty 
arising from the presence of an indefinite metric. The field schemes considered 
can be applied to a construction of quark models or for the classification of 
hadrons by non-compact groups (see [79]), and possibly, for studying sub- 
quark matter. 

The calculated density matrices (matrix-dyads) for fields and the method 
of computing the traces of 16— dimensional Petiau-Duffin-Kemmer matrix 
products allow us to make evaluations of different physical quantities in a 
covariant manner. 



A method of computing the traces of 16— dimensional Petiau-Duffin- 
Kemmer matrix products will be considered [80]. 

In order to compare the 16— component model of vector fields (including 
scalar states) with the Proca and Petiau-Duffin-Kemmer theories we consider 
here the density matrices in the form of matrix-dyads (40) for pure spin 
states. Taking into account Eqs. (16) it is possible to have the simpler 
expressions by using equalities 



where pW = p ^), p W = p^W, p (o) = p<°) = p^f\ The 



relationships (137) are obtained by using Eqs. (16), (24) and the expression 
for the squared spin operator: 



APPENDIX A 




(137) 




(138) 
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Taking into account relations (137), the projection matrix-dyads (40) are 
transformed into 

A W = -jLipW (ipW - em) a™ (a™ + s p ) = • (139) 

A W = -±- 2 ip& (ip& - em) ap (a® + Sp ) = • ^ {T) , (140) 

A?> = -L^ 1 ' (<p (1) - em) (l - < )2 ) = *?> • *?\ (141) 
A® = ^P (T) (V T) - em) (l - <) 2 ) = V® ■ ¥j\ (142) 



A (°) = ^tp {0) (V 0) - em) (l - af 2 ) 

= -Up (0) fe (0) - em) = ^ (0) • ^ {0) , 
2m 2 v ' 

1(0) = _L^(o) ^(o) _ £m ) ^ _ ^(0 



l_ ip (o)( ip (o) 



(143) 



(144) 



2m 2 

where we also used here the equalities: 

p {0) (i - < )2 ) = p (0 \ P (1 V P = = pW^^ 1 )^), 

I P I 

and so on. 

The matrix-dyads (139), (141) correspond to the vector state, (140), (142) 
- to the pseudovector state, (143) - to the scalar state, and (144) - to the 
pseudoscalar state. Eventually Eqs. (139)-(144) coincide with solutions of 
10— component (for vector and pseudovector states) and 5— component (for 
scalar and pseudoscalar states) free Petiau-Duffin-Kemmer equations. 

It is known that cross-sections for the scattering processes are summed 
to evaluate the transition probabilities for a particle going from the initial to 
the final states. These probabilities are proportional to the squared module's 
of the matrix elements which can be written as (see for example [52, 53]): 

| M | 2 =e 2 tr{QIIiQII 2 }, (145) 

where Q is the vertex operator, Q = i]Q + r] (77 = I^I^ is the Hermitian- 
izing matrix, Q + is the Hermite conjugated operator), the matrix-dyads ITi, 
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n 2 correspond to initial and final states, respectively. Therefore we need the 
traces of 16— dimensional Petiau-Duffin-Kemmer matrix products to calcu- 
late some processes with the presence of vector and scalar fields. 

It is easy to verify that the property of traces of the 16 x 16— Petiau- 
Duffin-Kemmer matrices (3^ (16) 

tr{/V/V--/V} = tr {(PP^P) (PP,A 3 P) {PP^-A^P)} 

(146) 

+tr {(P^ 2 P) (PM.P) - (P^n-AnP)} 

is valid, where P = e°'° + e°'° + (l/2)e^'^ is the projection operator. The 
analogous identity was derived in [48] for the 5 x 5— and 10 x 10— Petiau- 
Duffin-Kemmer matrices. The trace of the odd-numbered matrices (3^ is 
equal to zero. Eq. (146) is valid for any Petiau-Dumn-Kemmer matrix given 
by (16). Using the properties of the entire matrix algebra e A,B we find 

P^P^P = £ [H <° + \e^ P J^\ 

PPi +) Pi +) P = <W e °'° + £ M ' W , 
PPPP^P = 5, u e ' + \ex^e Xvaa s^™\ (147) 
Ppf^P = P^^P = e^\ 
Pti +) Pi 0) P = PPPPPP = ^ M '°, 

pW )p = pppffip = 

The relations (146), (147) with the help of the equalities tr|e A ' B | = 5a,b, 

= ^ p 5 ua - S^Svp, e AB e c ' D = 5 B c£ A,D , ( £AB ) CD = ^ac^bd allow us 
to find traces of any Petiau-Duffin-Kemmer matrices. 

APPENDIX B 
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We analize here the Lorentz transformations in the framework of the 
quaternion algebra. The quaternion algebra is defined by four basis elements 
e M = (ek,e±) (see, for example [81]) with the multiplication properties: 

2 1 2 2 2 1 

e 4 = 1, e 1 = e 2 = e 3 = -1, eie 2 = e 3 , 

e 2 e 1 = -e 3 , e 2 e 3 = e u e 3 e 2 = -e u (148) 

e 3 ei = e 2 , e x e z = -e 2 , e 4 e m = e m e 4 = e m , 

where m — 1, 2, 3 and e 4 = 1 is the unit element. 
The complex quaternion (or biquaternion) q is 

9 = 9^ = q m e m + q 4 e 4 , (149) 

where the q^ are complex numbers (q^ =Re q^ + ilm q^, i 2 = —1). Using 
the laws of multiplication (148), we find that the product of two arbitrary 
quaternions, q, q', is defined by: 

qq' = Mi - qmq'm) e 4 + (?4?m + 949m + ^mnkq n q'k) ^m- (150) 

It is convenient to represent the arbitrary quaternion as q = q 4 + q (so 
5 4 e 4 — > q 4 , q m e m — > q), where q 4 and q are the scalar and vector parts of the 
quaternion, respectively. With the help of this notation, Eq. (150) can be 
rewritten as 

qq' = q^q'i - q m q' m + qW + 94q' + [q, q'] • ( 15 i) 

Thus the scalar (q, q') = q m q' m i an< l vector [q, q'] products are parts of 
the quaternion multiplication. It is easy to verify the combined law for three 
quaternions: 

(9102)93 = 51(3253)- (152) 

The operation of quaternion conjugation (hyperconjugation) denotes the 
transition to 

q = q±e 4 - q m e m = q 4 - q, (153) 

so that the equalities 

91 + 92 = 9i + 9 2? 9i92 = 9 2 9i (154) 

are valid for two arbitrary quaternions q\ and q 2 . The modulus of the 
quaternion | q | is defined by 

\ q \= V(fi = \fqj- (155) 
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This formula allows us to divide one quaternions by another, and thus the 
quaternion algebra includes this division. For example, the simple equations 
and their solutions are given by 

q 2 x = q u x L - 



I 2 ' 

(156) 

xq 2 = qi, x R 



I Q2 



1 <?2 I 2 ' 

Quaternions are a generalization of the complex numbers and we can 
consider quaternions as a doubling of the complex numbers. They are con- 
venient for investigating the symmetry of fields and relativistic kinematics. 
In particular, the finite transformations of the Lorentz eigengroup are given 
by [81]: 

x' = LxL*, (157) 

where x = X4 + x is the quaternion of the coordinates (x^ = it, t is the 
time, x m are the spatial coordinates), L is the quaternion of the Lorentz 
group with the constraint LL = 1, L = L\ — L* and * means the complex 
conjugation. The biquaternion L with the constraint LL = 1 is defined by 
six independent parameters which characterize the Lorentz transformations. 
The squared four- vector of coordinates, x 2 , is invariant under the transfor- 
mations (157): 

x „ = x x = LxL L xL = xx = x^, (158) 

as L*L* = L*L* = 1. Eq. (158) shows that the 6— parameter transforma- 
tions (157) belong to the Lorentz group SO(3, 1). 

The complex Lorentz group SO (4, c) (see [82]) acts in complex space-time 
with the coordinates = (z m ,izo). The transformations of the complex 
Lorentz group are defined by 

z' = LzR, (159) 

where independent biquaternions L and R satisfy the requirement that 
LL = 1, RR = 1. As a result there are 12 independent parameters defining 
the SO(4:,c) group in which remains z^ invariant. Indeed, 

z* = z'z' = LzRRzL = zz = z\. (160) 

In the case of the ordinary Lorentz group we should put R = L* . 
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It should be noted that quaternion algebra can be realized using the Pauli 
2 x 2— matrices r k (k = 1, 2, 3). Setting e 4 = J 2 , e k = ir k and using the 
properties of Pauli's matrices: 

TmTn ^mnk^k ~\~ $mm 

(161) 

T%Tk + T k Ti = 2S ik , 

we come to the quaternion algebra (148). 
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